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Abstract 

An  industrial  burner  operating  in  the  MILD  combustion  regime  through  internal  recirculation  of  exhaust  gases 
has  been  characterized  numerically.  To  develop  a  self-sufficient  numerical  model  of  the  burner,  two  subroutines  are 
coupled  to  the  CFD  solver  to  model  the  air  preheater  section  and  heat  losses  from  the  burner  through  radiation.  The 
resulting  model  is  validated  against  experimental  data  on  species  concentration  and  temperature.  A  3 -dimensional 
CFD  model  of  the  burner  is  compared  to  an  axisymmetric  model,  which  allows  considerable  computational  saving, 
but  neglects  some  important  burner  features  such  as  the  presence  of  recirculation  windows.  Errors  associated  with 
the  axisymmetric  model  are  evaluated  and  discussed,  as  well  as  possible  simplified  procedures  for  engineering 
purposes.  Modifications  of  the  burner  geometry  are  investigated  numerically  and  suggested  in  order  to  enhance 
its  performances.  Such  modifications  are  aimed  at  improving  exhaust  gases  recirculation  which  is  driven  by  the 
inlet  air  jet  momentum.  The  burner  is  found  to  produce  only  30  ppmv  of  NO  when  operating  in  MILD  combustion 
mode.  For  the  same  air  preheating  the  NO  emissions  would  be  of  approximately  1000  ppmv  in  flame  combustion 
mode.  It  is  also  shown  that  the  burner  ensures  more  homogeneous  temperature  distribution  in  the  outer  surfaces 
with  respect  to  flame  operation,  and  this  is  attractive  for  burners  used  in  furnaces  devoted  to  materials’  thermal 
treatment  processes.  The  effect  of  air  excess  on  the  combustion  regime  is  also  discussed. 

©  2007  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 
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1.  Introduction 

The  improvement  of  combustion  efficiency  to 
reduce  fossil  fuel  consumption  and  carbon  diox¬ 
ide  emission  is  a  key  issue  in  combustion  research. 
Therefore,  the  achievement  of  high-energy-efficiency 
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processes  through  advanced  heat  recovery  systems  is 
desirable  even  if  the  increase  in  process  temperature 
arising  from  air  preheating  generally  leads  to  very 
large  NOx  emissions.  However,  NO*  emissions  are 
regulated  in  many  countries  with  increasingly  strin¬ 
gent  laws  because  of  their  adverse  impact  on  the 
environment  (climate  change,  acid  rains,  photochem¬ 
ical  smog).  Thus,  the  development  of  a  combustion 
technology  able  to  accomplish  large  energy  savings 
with  very  low  pollutant  emissions  has  become  a  pri¬ 
mary  concern. 
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area .  ur 

air  inlet  cross-sectional  area .  m 

Damkohler  number .  - 

emissivity .  - 

air  excess .  - 

view  factor .  - 

insulation  layer  thermal 

conductivity .  W  m-  1  K- 1 

kinetic  constant  of  prompt  NO _  s-1 

recirculation  degree  of  burnt  gases  into 

the  reaction  region .  % 

parameter  in  Eq.  (12) .  - 

kinetic  constant  of  thermal 

NO .  m1,5  K0,5  mol-0,5  s-0,5 

kinetic  constant  integrated  over  a  PDF 
of  temperature,  arbitrary  units 


emitter  length .  m 

mass  flow  rate .  kgs-1 

probability  density  function .  - 

irradiative  heat  flux .  W 

burner  load .  W 


heat  flux  exchanged  in  the  air 

preheater .  W 

radial  coordinate .  m 

radius .  m 

temperature .  K 

_ Q  _  1 

reaction  rate .  kg  m  J  s 

molecular  weight .  kg  mol- 1 


x  axial  coordinate .  m 

X  molar  fraction .  - 

[]  molar  concentration .  molm-3 

Greek  symbols 

tj  efficiency .  - 

p  density .  kgm-3 


o  Boltzmann  constant,  o  —  5.67  x 

10-8  W m-2  K-4 

rc  chemical  timescale .  s 

Tt  turbulent  timescale .  s 

Subscripts 

1  radiant  tube 

2  Inconel  shield 

2!  internal  insulation  layer  of  the  Inconel 

shield 

2!'  external  insulation  layer  of  the  Inconel 

shield 

3  water  heat  exchanger 

A  air 

EG  exhaust  gases 

F  fuel 

max  maximum 

mix  mixture 

rad  radiant  tube 

reaction  reaction  region 
rec  heat  recovery  section 

A k  parameter  in  Eq.  (12) 


In  this  framework,  MILD  combustion  appears  in¬ 
teresting,  as  it  may  ensure  high  combustion  efficien¬ 
cies  with  low  pollutant  emissions  [1,2]. 

Such  a  combustion  regime  needs  the  reactants  to 
be  preheated  above  the  self-ignition  temperature  and 
enough  inert  combustion  products  to  be  entrained  in 
the  reaction  region.  The  former  requirement  ensures 
high  thermal  efficiency,  whereas  the  latter  allows  di¬ 
luting  the  flame  and  reducing  the  final  temperature 
well  below  the  adiabatic  flame  temperature.  As  a  re¬ 
sult,  a  flame  front  is  no  longer  identifiable,  so  that 
MILD  combustion  is  often  denoted  as  flameless  com¬ 
bustion.  Moreover,  it  has  been  observed  that  igni¬ 
tion  and  extinction  phenomena  do  not  occur  in  MILD 
combustion  because  of  the  small  temperature  differ¬ 
ence  between  burnt  and  unburnt  gases. 

Among  the  advantages  associated  with  this  com¬ 
bustion  technology,  one  is  that  flame  stabilization  oc¬ 
curs  naturally  as  the  reactants’  temperature  exceeds 
the  self-ignition  temperature.  Therefore  a  large  de¬ 
gree  of  freedom  in  the  choice  of  the  fluid  dynamical 
configuration  of  the  combustion  chamber  is  allowed. 
Moreover,  the  temperature  field  homogeneity  and  re¬ 


duced  gradients  allow  better  control  of  maximum 
temperatures  with  beneficial  effects  on  materials. 

Importantly,  NO*  emissions  are  greately  reduced 
because  of  the  limited  temperature  increase.  Soot  for¬ 
mation  is  also  suppressed,  because  of  the  lean  con¬ 
ditions  in  the  combustion  chamber,  due  to  the  large 
dilution  levels.  In  addition,  the  large  CO2  concentra¬ 
tion  due  to  the  recirculation  of  combustion  products 
has  a  beneficial  effect  of  soot  suppression  [2]. 

Therefore,  MILD  combustion  poses  itself  as  a 
technology  combining  high  efficiencies,  because  of 
the  strong  preheating,  with  low  pollutant  emissions. 

All  these  aspects  make  MILD  combustion  worthy 
of  further  investigations  and  attention,  even  though,  to 
date,  several  studies  have  been  devoted  to  understand¬ 
ing  its  operational  conditions  [3]  as  well  as  its  mecha¬ 
nisms  and  critical  parameters  [4].  An  extensive  review 
on  MILD  combustion  features  considering  physical, 
chemical,  and  thermodynamic  aspects  has  been  pro¬ 
vided  by  Cavaliere  and  de  Joannon  [2]. 

From  a  technological  point  of  view,  the  first  re¬ 
quirement  for  MILD  combustion,  reactant  temper¬ 
ature  above  the  self-ignition  temperature,  may  be 
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achieved  by  preheating  the  fuel,  the  oxidizer,  or  both. 
The  second  requirement,  large  entrainment  of  inert 
species  in  the  reaction  region,  may  be  achieved  in  dif¬ 
ferent  ways  by  either  internal  or  external  recirculation 
of  exhaust  gases. 

MILD  combustion  burners  appear  well  suited  for 
those  industrial  processes  that  require  a  high  and  ho¬ 
mogeneous  temperature  distribution  within  the  com¬ 
bustion  chamber  (e.g.,  in  the  glass  and  ceramic  indus¬ 
try,  in  steel  thermal  treatments).  For  these  processes, 
energy  recovery  represents  a  primary  issue  to  ensure 
acceptable  energy  efficiency;  consequently,  the  use  of 
MILD  combustion  seems  to  be  particularly  benefi¬ 
cial  [1].  Energy  saving  can  be  achieved  with  either 
recuperative  or  regenerative  burners  according  to  the 
heat  exchanger  area  and  the  preheating  rate  required. 
These  burners  are  usually  designed  to  operate  both 
in  flame  and  in  MILD  mode  by  varying  the  reac¬ 
tants’  feeding  mode.  The  combination  of  recupera¬ 
tive  MILD  combustion  burners  with  radiant  tubes  is 
a  widely  applied  solution  in  many  thermal  treatments 
of  material  surfaces,  to  avoid  any  contact  or  contam¬ 
ination  of  the  flue  gases  with  the  stock  surface  to  be 
treated  [5]. 

Choi  and  Katsuki  [6]  investigated  the  feasibility 
of  flameless  oxidation  in  industrial  glass  furnaces. 
They  found  that  the  combustion  process  was  sus¬ 
tained  even  with  low-calorific-value  fuels  and  low 
oxygen  concentrations  if  the  combustion  air  was  pre¬ 
heated  above  the  fuel  self-ignition  temperature.  Re¬ 
sults  also  showed  that  NO*  formation  was  controlled 
by  the  mixing  process  between  fuel  and  the  preheated 
air. 

Flamme  [7]  investigated  the  application  of  MILD 
combustion  to  glass  melting  furnaces  operating  with  a 
process  temperature  around  1600  °C  and  air  preheat¬ 
ing  temperatures  up  to  1350  °C.  Experimental  results 
on  a  300-kW  furnace  showed  that  the  replacement  of 
conventional  burners  with  recuperative  MILD  com¬ 
bustion  burners  lead  to  dramatic  reduction  of  NO* 
emission,  from  1500  ppmv  down  to  values  safely  be¬ 
low  100  ppmv. 

The  application  of  MILD  combustion  in  different 
fields  than  thermal  treatment  processes  is  also  very 
attractive.  The  advantages  of  a  clean  and  quiet  com¬ 
bustion  process  could  be  exploited  in  several  fields 
of  application,  including  power  generation,  micro  co¬ 
generation,  and  low-temperature  applications. 

Flamme  [8]  showed  the  applicability  of  recupera¬ 
tive  MILD  combustion  burners  to  gas  turbines  in  over¬ 
coming  the  oscillation  problems  typical  of  lean  pre¬ 
mixed  combustion  occurring  in  gas  turbines.  Recently 
Wang  et  al.  [9]  performed  a  technical,  environmental, 
and  economic  analysis  of  NO*  reduction  technolo¬ 
gies  at  a  gas  turbine  power  station.  Results  showed 
that  the  use  of  recuperative  MILD  combustion  burn¬ 


ers  could  offer  an  effective  method  for  reducing  NO* 
emissions,  the  technology  being  much  cheaper  than 
SCR  and  requiring  only  a  slight  increase  in  capital 
cost  and  electricity  selling  price  with  respect  to  the 
simple  cycle  gas  turbine  power  plants. 

Also,  the  homogeneous  charge  compression  igni¬ 
tion  (HCCI)  engines  are  based  on  a  concept  that  may 
be  included  in  MILD  combustion  [2].  HCCI  engines 
achieve  a  homogeneous  charge  in  the  combustion 
chamber  at  such  a  high  compression  ratio  that  self¬ 
ignition  occurs  homogeneously  in  the  entire  cham¬ 
ber.  Under  this  condition  the  mixture  temperature  is 
higher  than  the  self-ignition  temperature;  however, 
the  maximum  temperature  is  kept  low  by  means  of 
a  superlean  mixture  and/or  inert  dilution.  With  such 
condition,  very  little  NO*  is  produced  and  relatively 
high-efficiency  cycles  can  be  designed. 

Given  the  industrial  interest,  recently  increasing 
attention  has  been  paid  not  only  to  the  experimen¬ 
tal  characterization  of  MILD  combustion  burners,  but 
also  to  their  modeling  through  computational  fluid 
dynamics  (CFD).  In  particular,  CFD  may  help  in  op¬ 
timizing  burners’  performances  by  investigating  geo¬ 
metrical  details,  such  as  injection  nozzles’  configu¬ 
rations  and  internal  devices  for  flue  gas  recircula¬ 
tion.  Moreover,  data  on  lab-scale  burners  are  dif¬ 
ficult  to  transfer  to  industrial  burners  operating  in 
MILD  combustion,  as  the  usual  scaling  criteria  are 
no  longer  adequate.  Indeed,  the  CFD  support  is  de¬ 
sirable  [10]. 

Some  attempts  at  modeling  MILD  combustion  op¬ 
eration  have  been  proposed  and  some  criteria  and 
hints  have  been  suggested,  even  if  there  is  still  not 
a  unique  response. 

Actually,  despite  the  flame  homogeneity  and  smo¬ 
oth  gradients  occurring  in  this  combustion  mode, 
MILD  combustion  appears  more  difficult  to  model 
than  conventional  flames,  e.g.,  premixed  or  diffusion 
flames,  because  the  high  dilution  levels  and  the  rel¬ 
atively  low  temperatures  lower  the  chemical  reac¬ 
tion  rates,  making  them  comparable  to  turbulent  mix¬ 
ing  phenomena,  which  are  enhanced  by  recirculation. 
Therefore  the  turbulence-chemistry  interaction  treat¬ 
ment  becomes  a  crucial  point  in  the  modeling  proce¬ 
dure. 

Coelho  and  Peters  [11]  simulated  a  furnace  oper¬ 
ating  in  the  MILD  combustion  mode  by  applying  the 
flamelet  approach  to  describe  turbulence/chemistry 
interactions.  They  compared  predictions  with  mea¬ 
surements  supplied  by  Pies  sing  et  al.  [12]  and  Ozde- 
mir  and  Peters  [13].  Although  a  qualitative  agreement 
between  the  experimental  and  predicted  flow  field 
was  observed,  some  considerable  discrepancies  were 
present.  In  addition  the  residence  time  was  underes¬ 
timated  by  the  numerical  model,  and  the  flamelet  ap¬ 
proach  was  found  to  be  unable  of  describing  correctly 
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the  NO  formation.  Better  results  could  be  achieved  by 
means  of  the  unsteady  flamelets  [14]. 

Orsino  et  al.  [15]  investigated  the  performance 
of  three  combustion  models  (eddy  break-up,  eddy 
dissipation  concept  with  chemical  equilibrium,  and 
PFD/mixture  fraction)  for  predicting  the  combustion 
of  natural  gas  with  high-temperature  air  and  large 
quantities  of  flue  gas  and  found  that  all  models  failed 
in  the  vicinity  of  the  fuel  injection.  However  the  au¬ 
thors  claimed  that  they  could  obtain  excellent  NO  pre¬ 
dictions  through  the  steady  flamelet  library. 

Recently,  Christo  and  Dally  [16]  pointed  out  dif¬ 
ficulties  arising  when  modelling  MILD  combustion 
burners.  They  investigated  numerically  a  jet  in  hot 
coflow  burner  replicating  MILD  combustion  regime. 
The  burner  was  characterised  experimentally  by  Dally 
et  al.  [17].  The  authors  used  different  turbulence, 
combustion,  and  kinetic  models,  comparing  their  per¬ 
formance.  In  particular  they  showed  that  mixture 
fraction/probability  density  function  and  flamelet  ap¬ 
proaches  perform  poorly  for  the  MILD  combustion 
regime.  The  eddy  dissipation  model  also  gave  un¬ 
satisfactory  results.  Better  predictions  were  achieved 
through  combustion  models  considering  both  chem¬ 
istry  and  turbulence  effects,  such  as  the  eddy  dissi¬ 
pation  concept  of  Magnussen  and  Hjertager  [18].  In 
addition,  the  aforementioned  authors  observed  that 
molecular  diffusion  may  play  an  important  role  in 
MILD  combustion:  for  instance,  the  predicted  tem¬ 
peratures  were  higher  than  the  measured  ones  by  ap¬ 
proximately  10%  if  differential  diffusion  effects  were 
neglected.  It  is  worth  noting  that  the  general  view  in 
turbulent  flames  modelling  is  to  neglect  the  molecular 
diffusion  term  and  to  take  into  account  only  turbulent 
diffusion. 

In  the  present  work,  a  numerical  model  of  an  in¬ 
dustrial  burner  operating  in  MILD  combustion  is  vali¬ 
dated  against  experimental  data.  Local  measurements 
across  the  combustion  chamber  could  not  be  per¬ 
formed  because  the  industrial  burner  does  not  allow 
any  access  to  its  interior;  however,  efforts  were  made 
to  characterize  the  burner  behavior  through  as  many 
measurements  as  possible  (for  instance,  by  measuring 
temperatures  along  the  outer  surfaces  of  the  burner). 
The  proposed  numerical  model  contains  two  subrou¬ 
tines  able  of  representing  the  recuperative  section  be¬ 
havior  and  the  heat  exchanges  between  the  burner  and 
the  surroundings.  These  subroutines  were  necessary 
in  order  to  avoid  any  assumption  in  the  burner  model¬ 
ing.  In  addition,  two  models  of  the  burner  have  been 
investigated:  a  3-dimensional  model,  which  replicates 
the  real  burner  features,  and  an  axisymmetric  model. 
The  axisymmetric  model  has  been  investigated  as 
this  model  is  commonly  adopted  [19,20]  because  of 
the  low  computational  cost.  However  the  axisymmet¬ 
ric  model  neglects  some  burner  features  that  could 


Fig.  1.  (a)  Longitudinal  section  of  the  burner  and  (b)  3-di¬ 
mensional  view  of  the  combustion  chamber. 


considerably  affect  the  fluid  dynamical  behavior  of 
the  burner  and,  subsequently,  predictions.  Ultimately, 
burner  geometry  modifications  and  effects  of  some 
operational  conditions  (e.g.,  air  excess)  on  the  burner 
behavior  are  numerically  analyzed. 


2.  The  burner 

The  burner  is  a  recuperative  MILD  combustion 
burner  with  a  nominal  power  of  13  kW  and  belongs  to 
ENEL  Ricerca.  A  longitudinal  section  of  the  burner 
is  reported  in  Fig.  la.  The  combustion  chamber  is 
delimited  by  a  radiant  tube  closed  at  the  upper  end. 
Indeed,  the  burner  is  suited  for  all  applications  where 
the  combustion  environment  has  to  be  kept  separated 
from  the  media  to  be  heated  (e.g.,  furnaces  for  steel 
annealing  processes,  glass  making).  The  combustion 
chamber  is  cylindrical  with  a  radius  of  0.045  m  and 
a  length  of  0.58  m.  It  operates  with  an  internal  re¬ 
circulation  of  exhaust  gases  which  is  promoted  by  a 
0.41-m-long  flame  tube  with  a  radius  of  0.02  m  posi¬ 
tioned  inside  the  burner.  The  flame  tube  is  equipped 
with  three  windows  in  the  lower  part  which  allow  the 
recirculation  of  exhaust  gases  in  the  reaction  region. 
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A  3 -dimensional  view  of  the  combustion  chamber 
with  the  recirculation  windows  is  depicted  in  Fig.  lb. 

Outside  the  combustion  chamber,  and  coaxial  to 
the  radiant  tube,  an  Inconel  shield  and  a  water  heat 
exchanger  were  placed  in  the  experimental  appara¬ 
tus  to  replicate  heat  losses  toward  the  surrounding  of 
real  burner  operations.  Two  insulation  layers  covered 
the  inner  and  outer  surfaces  of  the  Inconel  shield  (see 
Fig.  la). 

Importantly,  the  burner  is  self-recuperative,  which 
means  that  the  inlet  air  is  preheated  with  exhaust 
gases  by  means  of  a  finned- surface  heat  exchanger 
(see  the  lower  part  of  Fig.  la). 

The  burner  can  work  in  either  flame  or  MILD 
combustion  mode.  Initially,  the  combustion  chamber 
should  be  heated  up  with  a  flame  combustion  mode, 
as  MILD  combustion  requires  temperatures  higher 
than  the  self-ignition  temperature  of  the  fuel.  Once 
the  high  temperatures  are  reached,  the  burner  switches 
from  flame  to  MILD  combustion  mode. 

Fuel  and  air  are  fed  into  the  combustion  cham¬ 
ber  through  separated  jets.  In  particular,  the  fuel  is 
injected  on  the  burner  axis,  whereas  the  air  is  fed 
from  an  outer  annulus  in  a  co-flow  configuration.  In 
this  manner  the  fresh  air  can  entrain  large  amounts  of 
burnt  gases  before  reacting  with  the  fuel. 

Measurements  were  made  of  the  composition  of 
exhaust  gases  as  well  as  temperatures  along  the  radi¬ 
ant  tube  and  the  Inconel  shield,  both  of  them  equipped 
with  series  of  thermocouples. 

Experimental  runs  are  listed  in  Table  1  (see  runs 
1-7).  Different  burner  loads  were  investigated  by 
varying  the  fuel  flow  rate  and  air  excess. 

3.  Numerical  model 

The  numerical  simulations  were  performed  with 
the  commercial  code  CFX  5.7  by  Ansys  Inc. 


It  can  be  observed  from  Table  1  that  runs  1-7 
were  investigated  both  numerically  and  experimen¬ 
tally,  whereas  runs  8-14  were  analyzed  only  numer¬ 
ically.  These  latter  runs  were  partly  aimed  at  investi¬ 
gating  burner  geometry  modifications  (runs  8-12)  and 
the  influence  of  air  excess  (runs  13-14)  on  the  burner 
behavior. 

In  particular,  the  air  inlet  cross-sectional  area  was 
varied  from  339  down  to  41  mm  ,  in  order  to  change 
the  recirculation  degree  of  exhaust  gases,  such  recir¬ 
culation  being  driven  by  the  momentum  of  the  air 
inlet  jet.  Actually,  the  real  burner  operates  with  just 
one  air  inlet  cross-sectional  area,  88  mirr.  Values  of 
A  air  ?  in  f°r  each  run  are  listed  in  the  fifth  column  of 
Table  1. 

The  air  excess  was  varied  from  40  to  46%,  which 
characterizes  the  experimental  campaign  and  ensures 
an  oxygen  concentration  of  about  6%  in  the  exhaust 
gases,  down  to  16%,  corresponding  to  an  oxygen  con¬ 
centration  of  about  2.6%  in  the  exhaust  gases. 

3.1.  Computational  domain  and  grid 

The  computational  domain  consisted  of  three  sub- 
domains:  a  fluid  domain  representing  the  real  com¬ 
bustion  chamber  and  two  solid  domains  represent¬ 
ing  the  flame  and  the  radiant  tubes.  It  was  chosen  to 
simulate  those  two  solid  subdomains  to  take  into  ac¬ 
count  heat  exchanges,  due  mainly  to  radiation  effects 
between  the  reacting  gas  and  the  solid  boundaries. 
Besides,  the  radiant  tube  modeling  was  done  on  pur¬ 
pose  to  allow  a  comparison  between  predictions  and 
experimental  data  obtained  from  the  thermocouples 
mounted  along  the  radiant  tube. 

Because  of  the  presence  of  three  recirculation 
windows,  a  120°  angular  sector  was  modeled.  The 
grid  contained  350,000  cells  and  consisted  of  parts 
meshed  with  tetrahedra  (near  burner  zone)  and  parts 


Table  1 


Experimental  and  numerical  runs  investigated 


Run 

Gin  (kW) 

mF  (kg/s) 

^air  (%) 

^airjn  (mm2) 

£r  (%) 

C/O reaction  (  ) 

Ac>2,  reaction  (  ) 

Exp. 

CFD 

CFD  model 

1 

8.23 

1.72F-04 

60 

88 

130 

0.105 

0.126 

X 

X 

3D 

2 

9.13 

1.90F-04 

44 

88 

138 

0.122 

0.117 

X 

X 

3D 

3 

9.31 

2.38F-04 

40 

88 

137 

0.128 

0.115 

X 

X 

3D 

4 

9.9 

2.55F-04 

40 

88 

134 

0.129 

0.115 

X 

X 

3D 

5 

10.53 

2.70F-04 

44 

88 

135 

0.120 

0.118 

X 

X 

3D 

6 

10.7 

2.73F-04 

40 

88 

134 

0.129 

0.115 

X 

X 

3D 

7 

10.42 

2.67F-04 

46 

88 

136 

0.120 

0.119 

X 

X 

2D,  3D 

8 

10.42 

2.67F-04 

46 

339 

23 

0.160 

0.171 

X 

2D,  3D 

9 

10.42 

2.67F-04 

46 

201 

59 

0.144 

0.146 

X 

2D,  3D 

10 

10.42 

2.67F-04 

46 

64 

173 

0.111 

0.111 

X 

2D,  3D 

11 

10.42 

2.67F-04 

46 

46 

200 

0.105 

0.106 

X 

2D,  3D 

12 

10.42 

2.67F-04 

46 

41 

235 

0.098 

0.102 

X 

2D,  3D 

13 

10.42 

2.67F-04 

16 

88 

128 

0.183 

0.099 

X 

3D 

14 

10.42 

2.67F-04 

16 

41 

223 

0.165 

0.078 

X 

3D 
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Fig.  2.  (a)  3 -Dimensional  and  (b)  axisymmetric  grids. 


meshed  with  hexahedra  (reaction  and  reverse  flow 
zone,  flame  tube  and  radiant  tube).  Specifically  there 
were  143,000  tetrahedra  and  207,000  hexahedra.  The 
choice  of  a  hybrid  grid  was  aimed  at  reducing  the 
number  of  cells  required  to  discretize  the  computa¬ 
tional  domain  accurately.  Tetrahedra  were  required  to 
properly  characterize  some  burners’  features,  e.g.,  re¬ 
circulation  windows;  however,  a  grid  constituted  only 
of  tetrahedra  would  have  been  characterized  by  a  huge 
number  of  cells.  The  computational  grid  is  shown  in 
Fig.  2a. 

Some  simulations  were  also  performed  with  the 
axial  symmetry  hypothesis,  thus  neglecting  the  pres¬ 
ence  of  the  recirculation  windows.  In  other  words,  in 
these  simulations  the  recirculation  area  corresponded 
to  the  whole  annulus  comprised  between  the  lower 
part  of  the  burner  and  the  flame  tube;  consequently 
the  recirculation  area  is  larger  than  that  considered 
with  the  3D  model.  This  was  done  on  purpose  to  eval¬ 
uate  the  effect  of  recirculation  windows;  nevertheless, 
simulations  available  in  the  literature  of  recuperative 
MILD  combustion  burners  are  usually  2D  simulations 
[19,20].  The  axisymmetric  model  is  shown  in  Fig.  2b. 
The  structured  grid  consisted  of  120,000  hexahedra, 
arranged  in  four  layers  describing  an  angular  sector 
of  5°. 

A  grid  independence  study  was  carried  out  on  both 
the  3D  and  the  axisymmetric  models.  The  effect  of 
mesh  refinement  and  coarsening  was  investigated  by 
doubling  and  halving  the  number  of  elements  of  the 


actual  mesh  used  in  the  computations  and  comparing 
the  radial  profiles  of  temperature  and  major  species  at 
different  axial  locations  along  the  burner. 

3.2.  Boundary  conditions 


As  far  as  boundary  conditions  are  concerned,  fu¬ 
els  were  fed  at  room  temperature  (298  K).  Conversely, 
the  air  inlet  temperature  was  evaluated  as  a  function  of 
exhaust  gas  temperature,  on  the  basis  of  experimen¬ 
tal  data  on  the  air  preheater  available  from  the  burner 
supplier.  A  model  was  developed  from  experimental 
data  fitting  and  implemented  through  a  subroutine  in 
the  code;  as  a  result,  the  air  inlet  temperature  was  ad¬ 
justed  on  the  basis  of  the  exhaust  gas  temperature. 

As  mentioned  earlier,  special  attention  was  paid  to 
the  boundary  conditions  of  the  radiant  tube,  as  some 
of  the  available  experimental  data  regarded  this  re¬ 
gion.  Heat  is  lost  from  the  radiant  tube  toward  the 
Inconel  shield  by  radiation,  and  subsequently  heat 
is  conducted  through  the  Inconel  shield  and  the  two 
insulation  layers  covering  it;  ultimately  heat  is  trans¬ 
ferred  from  the  shield  to  the  water  heat  exchanger 
through  radiation. 

A  subroutine  based  on  conduction  and  radiation 
between  coaxial  cylindrical  shields  was  written  on 
purpose  to  evaluate  the  heat  losses  from  the  radi¬ 
ant  tube.  The  irradiative  heat  flux  between  radiant 
tube  (1)  and  inner  insulation  layer  (2r)  is 


Qn'  =  en'°^\0\-Tl3 


(1) 


where  eyi'  —  — 


1 


e\ 


l 


^1  (  1  e2r 


■(— 21) 
V  / 


el  F\2r  R2'  v  e2' 

The  conductive  heat  transfer  within  the  insulation 

>/ _ A  o// 


layers  2  and  2  is  given  by 


27rLem/qns(7V  —  ^2") 

ln(/?2"/^2') 


neglecting  the  thermal  resistance  of  the  Inconel 
shield. 

The  irradiative  heat  flux  between  the  outer  insula¬ 
tion  layer  (2")  and  the  water  heat  exchanger  (3)  is 


2 2" 3  —  e2"2G  ^2"  {T2"  ~  ^3)  ’  (3) 

where 


1 


In  the  above  equations  T\  is  the  radiant  tube  tem¬ 
perature  and  73  is  the  water  heat  exchanger  temper¬ 
ature,  i.e.,  73  =  323  K,  whereas  T%  and  7^//  are  the 
temperatures  of  the  inner  and  outer  insulation  layers 
respectively,  e^ey,  ^2"»  and  e3  are  emissivities  of 
the  radiant  tube,  the  inner  and  outer  insulation  layers 
of  the  Inconel  shield,  and  the  water  heat  exchanger, 
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respectively.  Lem  is  the  emitter  length,  i.e.,  the  length 
of  the  coaxial  shield  radiating  toward  the  surround¬ 
ings,  whereas  k[ns  is  the  thermal  conductivity  of  the 
insulation  layers.  F\2f  and  Fy 3  are  the  view  factors 
between  coaxial  cylinders;  thus  their  values  were  set 
equal  to  1  (neglecting  radiation  losses  towards  the 
burner  axis). 

Such  a  value  is  strictly  correct  for  infinite  length 
cylinders;  however,  this  was  considered  an  adequate 
approximation  because  the  radiant  tube,  the  Inconel 
shield,  and  the  water  heat  exchanger  have  a  length 
much  greater  than  their  radius. 

At  steady  state  Q12'  —  Ql'l"  —  Q->  consequently 
the  following  expression  for  Q  is  obtained  by  express¬ 
ing  T2'  and  T2»  as  a  function  of  T\  and  Ty. 


Q  = 


1 


+ 


1 


oey2jA.\  g  e  2/ /  ^  -A  2 


x 


T\  -  T3  + 


T$ - 


Q 


0.25 


oe\2'A\ 


Q 

2  TC  TenAins 

In  R2ff/R2'  - 


—) 


Equation  (4)  made  necessary  an  iterative  procedure  to 
solve  the  heat  flux  Q . 

It  may  be  argued  that  a  discretization  of  the  radi¬ 
ant  tube  and  the  outer  cylinders  would  led  to  a  more 
accurate  evaluation  of  radiation  losses.  However,  the 
temperature  gradients  along  the  radiant  tube  were  not 
high  enough  to  justify  the  higher  computational  cost 
required  by  this  approach. 


3.3.  Physical  model 


influence  of  molecular  diffusion  on  flame  characteris¬ 
tics. 

The  discrete  transfer  radiation  model  was  used 
with  a  number  of  rays  equal  to  16.  A  coarsening  rate 
equal  to  24  was  applied  to  the  actual  grid  to  facilitate 
the  radiation  calculation.  Some  simulations  were  also 
made  with  the  Monte  Carlo  model,  but  as  results  dif¬ 
fered  by  less  than  5%,  the  former  model  was  preferred 
because  of  its  lower  computational  cost. 

The  radiation  properties  of  the  reacting  mixture 
were  taken  into  account  with  the  weighted- sum-of- 
gray-gases  (WSGG)  model  using  the  coefficients  pro¬ 
posed  by  Smith  et  al.  [22].  This  model  was  applied  by 
coupling  the  code  with  a  subroutine  written  on  pur¬ 
pose. 

As  far  as  the  kinetic  mechanism  is  concerned,  a 
simple  one- step  global  mechanism  was  used  for  the 
methane  oxidation  [23], 

^ch4 

CH4  +  202  — ►  C02  +  2H20,  (5) 

with 

^CH4  =  Wch48-3  X  105exp(-151°0K^ 

x  [CH4ra3[02] 

\mJ  s  / 

Thermal  NO  formation  was  modeled  with  a  finite 
rate  chemistry  combustion  (FRC)  model  and  a  sim¬ 
plified  one- step  kinetic  mechanism,  available  in  the 
code.  The  simplified  mechanism  is  obtained  from  the 
Zeldovich  scheme  by  assuming  a  steady  state  for  the 
N  radicals  and  relating  the  O  radical  concentration 
to  that  of  oxygen  by  means  of  the  dissociation  reac¬ 
tion  [24].  The  resulting  rate  is  expressed  as 


Favre- averaged  Navier-Stokes  equations  were 
solved  by  using  a  modified  version  of  the  k-s  tur¬ 
bulence  model.  The  first  constant  of  the  dissipation 
transport  equation  Ce\  was  set  equal  to  1.6  instead 
of  1.44  as  suggested  by  Morse  [21]  in  order  to  over¬ 
come  the  deficiency  of  the  standard  k-s  model  in 
predicting  round  jets  properly. 

The  combined  eddy  dissipation  model/finite  rate 
chemistry  (EDM/FRC)  was  chosen  as  combustion 
model.  Such  a  model  calculates  both  a  chemical  Ar¬ 
rhenius  rate  and  a  turbulent  mixing  rate  based  on  the 
Magnussen-Hjertager  equation  and  then  chooses  the 
lower  of  the  two  rates  to  be  inserted  in  the  species’ 
transport  equation.  The  EDM/FRC  approach  was  con¬ 
sidered  well  suited  for  modeling  MILD  combustion. 
As  mentioned  previously,  the  conditions  of  elevated 
and  uniform  temperature  distribution  and  low  oxygen 
concentration,  which  characterize  the  MILD  combus¬ 
tion  regime,  lead  to  slower  reaction  rates  with  respect 
to  conventional  combustion  regimes  and  enhance  the 


^NO, thermal  —  ^NO  ^thermal  [02J 


9 


where 


^thermal 

_  4.52  x  1015 

Cf 


exp  | 


69,466  K 


/  m3K  \1/2 

V  kmol  s2 


The  prompt  NO  formation  from  methane  was  mod¬ 
eled  using  a  similar  approach,  according  to  the  one- 
step  mechanism  proposed  by  De  Soete  [25], 


^NO, prompt  =  WNokprompt[02]  1//2[N2][F] 


where 


^prompt  —  1.2  x  10  exp  I 


30,215  K\  /I 
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For  both  thermal  and  prompt  NO  formations,  the  Ar¬ 
rhenius  equation  was  integrated  over  a  probability 
density  function  (PDF)  for  temperature  in  order  to 
take  into  account  turbulence  effects, 

max 

k(T)  =  - - h— -  [  k(T)P(T)dT ,  (8) 

1  max  1  min  J 

min 

where  rmax  —  Tm[n  is  the  temperature  range  of  in¬ 
terest  and  PDF(T)  is  the  probability  density  function 
that  is  computed  from  the  mean  temperature  and  the 
variance  of  the  temperature.  The  shape  of  the  PDF(T) 
is  presumed  to  be  a  two-moment  ^-function  as  appro¬ 
priate  for  combustion  calculations  [26].  Equation  (8) 
must  be  integrated  at  every  grid  node  and  at  every  it¬ 
eration.  The  limits  of  integration  are  determined  from 
the  minimum  and  maximum  values  of  temperature  in 
the  combustion  solution. 

The  implementation  of  detailed  NO  formation 
mechanisms  is  in  progress. 

Molecular  diffusion  effects  were  taken  into  ac¬ 
count  with  the  Slattery  and  Bird  [27]  equation  for 
binary  diffusion  coefficients.  The  molecular  species 
were  assumed  to  diffuse  in  a  mixture  manly  consti¬ 
tuted  by  nitrogen,  due  to  the  large  amount  of  inert 
gases  recirculating  in  the  reaction  zone.  This  approx¬ 
imation  was  found  to  be  consistent  with  results  ob¬ 
tained  from  a  more  rigorous  approach,  consisting  in 
the  evaluation  of  the  effective  mass  diffusivities  for 
the  species  diffusing  in  the  gas  mixture  by  applying 
the  Wilke  mixing  rule  [28]  to  the  binary  diffusion 
coefficients  obtained  from  the  Slattery  and  Bird  equa¬ 
tion  [27].  Thus,  the  former  approximation  was  pre¬ 
ferred  due  to  the  much  lower  computation  cost. 

However,  simulations  performed  with  and  without 
taking  into  account  molecular  diffusion  showed  that 
for  the  investigated  cases  the  effect  of  molecular  dif¬ 
fusion  is  negligible  with  respect  to  turbulent  diffusion. 

4.  Results 

4.1.  Assessment  of  the  numerical  model 

First,  CFD  simulations  were  validated  against  ex¬ 
perimental  data.  This  is  a  necessary  operation,  as  the 
CFD  model  of  the  burner  is  rather  complex,  contain¬ 
ing  two  subroutines  for  the  air  preheater  modeling  and 
the  calculation  of  heat  exchanges  between  the  burner 
and  surroundings  (see  Section  3.2). 

Fig.  3  shows  the  comparison  between  predicted 
and  measured  temperatures  along  the  radiant  tube  for 
two  different  experimental  runs  (runs  1  and  2  of  Ta¬ 
ble  1).  The  agreement  is  satisfactory,  as  the  temper¬ 
ature  predicted  profiles  fit  measured  data  well.  This 
confirms  that  the  representation  of  radiation  losses 


Fig.  3.  Comparison  between  predicted  and  measured  tem¬ 
peratures  along  the  radiant  tube  for  runs  1  and  2. 


Fig.  4.  Comparison  between  predicted  and  measured  NO 
emissions  normalized  to  3%  oxygen  in  the  flue  gases. 

through  a  model  based  on  radiation  between  coaxial 
cylindrical  shields  was  somewhat  appropriate,  despite 
the  approximations  behind  the  model. 

The  parity  plot  of  Fig.  4  shows  the  comparison 
between  experimental  and  predicted  NO  emissions. 
NO  concentrations  have  been  normalized  to  3%  of 
oxygen  in  the  exhaust  gases  in  order  to  compare  the 
different  runs,  these  being  characterized  by  different 
air  excesses.  It  can  be  observed  an  acceptable  agree¬ 
ment  between  predicted  and  measured  data,  although 
kinetics  of  NO  formation  were  rather  simple. 

4.2.  Effect  of  increasing  air  inlet  velocity 

The  flow  pattern  in  the  burner  confirms  the  mech¬ 
anism  that  promotes  the  exhaust  gases  recirculation. 
The  burnt  gases  are  entrained  by  the  fresh  air  jet  and 
deviated  from  their  path,  so  that  they  recirculate  to¬ 
ward  the  combustion  zone. 

The  degree  of  recirculation  can  be  evaluated  with 
the  recirculation  factor,  defined  as  =  fi?EG/(^B  + 
mpj,  where  m^Q  is  the  mass  flow  rate  of  the  exhaust 
gases  recirculating  into  the  reaction  zone,  k r  was 
found  to  increase  when  the  air  inlet  cross-sectional 
area  was  reduced,  varying  between  &r  =  23%  with 
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Fig.  5.  Temperature  distribution  for  (a)  Aajr  [n  =  339  mm2;  (b)  Aaa- =  201  mm2;  (c)  Aair  in  =  88  mm2;  (d)  Aair jn  =  64  mm2; 
(e)  Aajr  [n  =  46  mm2;  (f)  Aaa-  jn  =  41  mm2.  Burnerioaci  Qm  =  10.42  kW  (runs  7-12  of  Table  1). 


r\ 

A  air,  in  —  339  mm  and  k r  =  235%  with  — 

41  mm  ,  for  the  10.42-kW  burner  loads  (runs  7-12 
of  Table  1).  The  real  burner  operates  with  A^v  -m  = 
88  mm  ,  thus  with  a  recirculation  degree  of  k r  = 
136%.  Values  of  k r  for  each  run  are  listed  in  the  sixth 
column  of  Table  1 .  Moreover,  the  seventh  and  eighth 
columns  of  the  table  report  the  C/O  ratio  and  the  O2 


molar  fraction  in  the  reaction  region,  respectively.  It 
can  be  observed  that  the  C/O  ratio  is  rather  low,  but 
this  is  partly  imputed  to  the  relatively  large  air  ex¬ 
cesses  used  in  the  experimental  runs. 

Fig.  5  shows  the  temperature  distributions  in  the 
burner  for  different  recirculation  degrees.  A  more  ho¬ 
mogeneous  temperature  distribution  can  be  seen  in 
the  burner  with  increasing  recirculation  degree. 
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Fig.  6.  Axial  temperature  profiles  at  r  =  7  mm  for  different 
air  inlet  cross-sectional  areas.  Burner  load  Qm  =  10.42  kW 
(runs  7-12  of  Table  1). 

The  maximum  temperatures  observed  in  the  burner 
were  found  to  be  reduced  by  more  than  100  K  when 
the  air  inlet  cross-sectional  area  was  halved,  from 
Tmax  =  2080  K  for  &r  =  136%  down  to  Tin  ax  — 
1972  K  for  &r  =  235%.  Simulations  performed  with 
air  inlet  cross-sectional  areas  larger  than  the  real  one 
exhibited  high-temperature  peaks  of  2572  and  2351  K 
for  &r  =  23%  and  &r  =  59%,  respectively.  In  ad¬ 
dition,  these  predictions  indicated  the  presence  of  a 
flame  front;  thus  the  burner  was  operating  in  flame 
rather  than  in  MILD  combustion  conditions. 

Fig.  6  illustrates  axial  profiles  of  temperature  at  a 
radial  distance  r  —1  mm  and  for  different  air  inlet 
cross-sectional  areas.  The  radial  distance  r  =  1  mm 
was  chosen  as  it  corresponds  to  the  middle  of  the  inlet 
air  jet. 

It  can  be  noticed  that  the  temperature  increase 
is  by  about  850  K  for  the  real  burner  (Aa = 

r\ 

88  mmz),  whereas  it  is  reduced  down  to  750  K  for  the 
largest  recirculation  degree  (Aajr  =  41  mm~).  The 
methane  self-ignition  temperature  is  853  K,  therefore, 
according  to  the  criteria  by  Cavaliere  and  de  Joan- 
non  [2]  the  real  burner  operation  can  scarcely  be  re¬ 
garded  as  MILD  combustion  (“A  combustion  process 
is  named  MILD  when  the  inlet  temperature  of  the 
reactant  mixture  is  higher  than  mixture  self-ignition 
temperature  whereas  the  maximum  allowable  temper¬ 
ature  increase  with  respect  to  inlet  temperature  during 
combustion  is  lower  than  mixture  self-ignition  tem¬ 
perature  (in  Kelvin)”). 

For  smaller  recirculation  degree,  the  temperature 
increase  due  to  combustion  heat  release  is  drastic, 

r\ 

for  instance  being  of  1300  K  for  Aa =  339  mmz 
(£r  =  23%),  and  reasonably  the  burner  does  not  op¬ 
erate  in  MILD  combustion  mode. 

In  addition,  the  shape  of  temperature  profiles  ob¬ 
tained  for  &r  =  23  and  59%  differs  considerably  from 
that  obtained  with  larger  recirculation  degrees.  For 
large  recirculation  degrees  the  high  temperature  of 


recirculating  exhaust  gases  increases  the  mean  re¬ 
actant  temperature,  leading  to  temperature  profiles 
characterized  by  a  double  peak  (Fig.  6).  Indeed,  the 
first  temperature  peak,  located  at  x  =  100-150  mm, 
is  due  to  the  mixing  of  fresh  reactants  with  ex¬ 
haust  gases.  Conversely,  for  smaller  recirculation  de- 
grees  (Aair  =  201  and  339  mirr),  the  contribution 
of  exhaust  gases  to  the  reactants’  heat  capacity  is 
negligible,  so  that  the  first  temperature  peak  disap¬ 
pears. 

For  a  better  understanding  of  temperature  distrib¬ 
utions  in  the  MILD  combustion  burner,  it  is  useful  to 
analyze  radial  profiles  of  temperature  at  different  ax¬ 
ial  distances  along  the  burner,  as  reported  in  Fig.  7.  It 
can  be  observed  that  the  profiles  become  less  steeper 
with  increasing  the  recirculation  degree,  thus  con¬ 
firming  more  homogeneous  temperature  distributions. 
Again,  the  profiles’  shape  for  the  &r  =  23  and  59% 
runs,  corresponding  to  the  flame  combustion  mode, 
differs  strongly  from  that  of  the  MILD  combustion 
runs. 

4.3.  Effect  of  air  excess 

As  mentioned  in  the  Introduction,  the  experimen¬ 
tal  campaign  was  carried  out  with  a  relatively  high  air 
excess,  ensuring  an  oxygen  concentration  in  the  flue 
gases  of  6%.  For  that  reason,  simulations  were  per¬ 
formed  also  with  a  smaller  air  excess  of  16%,  leading 
to  an  oxygen  concentration  in  the  exhaust  gases  of  ap¬ 
proximately  2.6%  (see  runs  13  and  14). 

It  was  found  that  for  the  same  air  inlet  cross- 
sectional  area  Aa =  88  mm“,  the  maximum  tem¬ 
perature  observed  in  the  burner  increased  from  2080 
to  2265  K  when  the  air  excess  decreased  from  46 
to  16%.  This  may  be  easily  imputed  to  the  reduced 
thermal  capacity  associated  with  the  recirculating  ni¬ 
trogen.  In  addition,  NO  was  also  found  to  increase 
from  about  39  to  272  ppmv.  Even  if  the  oxygen  con¬ 
centration  in  the  reaction  zone  diminishes  with  de¬ 
creases  air  excess  (see  Table  1),  thus  damping  NO 
formation,  the  enhanced  temperatures  govern  the  NO 
formation  and  lead  to  a  strong  increase  of  NO  emis¬ 
sions.  Therefore,  with  an  air  excess  of  16%,  the  stan¬ 
dard  burner  configuration  is  not  sufficient  to  ensure  a 
completely  MILD  combustion  mode,  and  larger  recir¬ 
culation  degrees  would  be  desirable. 

With  the  smallest  air  inlet  cross-sectional  area, 

r\ 

A  air,  in  —  41  mm  ,  and  the  same  air  excess  of  16%, 
the  NO  emissions  were  found  to  be  reduced  down  to 
58  ppmv. 

4.4.  Turbulence-chemistry  interaction 

MILD  combustion  is  typically  characterized  by 
a  particular  kind  of  turbulence-chemistry  interac- 
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Fig.  7.  Radial  profiles  of  temperature  for  different  recirculation  degrees  under  different  axial  coordinates:  (a)  x  =  150;  (b)  250; 
(c)  350;  (d)  450  mm.  Burner  load  Qm  =  10.42  kW  (runs  7-12  of  Table  1). 


tion.  Specifically,  higher  turbulence  levels  with  re¬ 
spect  to  conventional  flames  occur  in  the  reaction 
region  because  of  the  strong  recirculation,  whereas 
slower  chemical  rates  are  observed  because  of  the 
large  dilution  of  reacting  species. 

A  useful  way  of  estimating  the  turbulence-chem¬ 
istry  interaction  is  provided  by  the  Damkohler  num¬ 
ber,  which  is  the  turbulence-to-chemical  timescale  ra¬ 
tio, 

tt 

Da  =  — .  (9) 

Fig.  8  shows  the  Damkohler  number  distribution  in 
the  burner  through  contours  for  different  recircula¬ 
tion  degree.  The  contours  have  been  overlapped  to  the 
temperature  distribution  in  the  burner  in  order  to  fa¬ 
cilitate  the  figure  comprehension.  Actually,  since  the 
chemical  timescale  is  based  on  one-step  kinetics,  re¬ 
sults  cannot  provide  a  quantitative  answer  but  give 
qualitative  evidence  of  the  effect  of  exhaust  gases  re¬ 
circulation  on  turbulence-chemistry  interactions. 

Fig.  8a  refers  to  the  largest  air  inlet  cross-sectional 
area,  for  which  the  burner  is  operating  in  flame  mode 
(£r  =  23%).  It  can  be  noticed  that  the  flame  front  is 
characterized  by  large  Da  values. 

Conversely,  in  Figs.  8b  and  8c,  that  Damkohler 
number  values  are  low  and  near  unity  even  in  the  re¬ 


action  region,  thus  proving  that  both  turbulence  and 
chemistry  control  the  reaction  rates.  This  indicates 
that  the  combustion  model  (EDM/FR)  used  in  the 
present  work  is  conceptually  consistent. 

Moreover,  it  can  be  observed  by  comparing  Fig.  8b 
and  Fig.  8c  that  regions  of  higher  Damkohler  num¬ 
bers  enlarge  with  increasing  the  recirculation  degree; 
this  is  due  to  the  larger  turbulence  and  reactants’  dilu¬ 
tion. 

4.5.  Burner  efficiency 

As  mentioned  in  the  introduction,  MILD  com¬ 
bustion  burners  are  very  interesting  from  an  indus¬ 
trial  point  of  view  and  have  been  applied  to  different 
processes.  Therefore  it  is  important  to  analyze  not 
only  their  “internal”  behavior  (e.g.,  temperature  and 
flow  fields)  but  also  their  interactions  with  the  sur¬ 
roundings  and  their  performance  in  terms  of  thermal 
efficiency.  In  particular,  recuperative  MILD  combus¬ 
tion  burners  enclosed  by  radiant  tubes  are  widely  used 
in  continuous  annealing  furnaces  for  electric  steel 
strip.  In  such  applications,  it  would  be  desirable  to 
achieve  smooth  temperature  profiles  along  the  radi¬ 
ant  tube  in  order  to  ensure  good  steel  quality  and 
safe  conditions  for  the  material  constituting  the  radi¬ 
ant  tube.  Therefore,  local  temperature  peaks  and  gra- 
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Fig.  8.  Temperature  (colored  map)  and  Damkohler  number  (contours)  distributions  for  (a)  Aajrjn  =  339  mm2;  (b)  Aajrjn 
88  mm2;  (c)  Aajrjn  =  46  mm2.  Burner  load  Q[n  =  10.42  kW  (runs  8,  7,  11  of  Table  1). 


dients  should  be  minimized  to  avoid  the  collapse  of 
the  radiant  tube  itself  due  to  localized  phenomenona, 
such  as  deformation  and  mechanical  stresses  arising 
from  nonuniform  thermal  expansion  [29] . 

Fig.  9  shows  temperature  profiles  along  the  radiant 
tube  for  a  burner  load  of  10.42  kW  and  different  val¬ 
ues  of  the  recirculation  degrees.  A  large  temperature 
interval,  greater  than  230  K,  is  observed  for  the  flame 
combustion  mode,  whereas  in  the  MILD  combustion 
regime  the  temperature  profiles  appear  more  uniform, 


covering  temperature  intervals  of  about  150  K.  There¬ 
fore  this  condition  is  more  suited  to  ensure  homoge¬ 
neous  conditions. 

Another  criterion  to  evaluate  the  burner  perfor¬ 
mance  is  based  on  the  thermal  efficiency.  In  particular 
two  different  efficiencies  can  be  calculated:  a  radia¬ 
tion  efficiency  based  on  heat  transfer  from  the  radiant 
tube  and  a  recovery  efficiency  characterizing  the  air 
preheater  performance.  The  former  efficiency  may  be 
calculated  as 


C.  Galletti  et  al.  /  Combustion  and  Flame  151  (2007)  649-664 


661 


Fig.  9.  Radiant  tube  temperature  as  a  function  of  the  air 
inlet  cross-sectional  area.  Burner  load  Qm  =  10.42  kW 
(runs  7-12  of  Table  1). 


Fig.  11.  Recirculation  degree  as  a  function  of  the  air  in¬ 
let  cross-sectional  area  predicted  with  the  axisymmetric  and 
3-dimensional  models.  Burner  load  Q\n  =  10.42  kW  (runs 
7-12  of  Table  1). 
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Fig.  10.  Preheater  and  radiant  rube  efficiencies  as  a  function 
of  the  recirculation  degree.  Burner  load  (9  in  =  10.42  kW 
(runs  7-12  of  Table  1). 
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whereas  the  recovery  efficiency  is 


(10) 


Oyqc  — 


Q  ICC 

£hT 


(11) 


Fig.  10  shows  the  two  efficiencies  as  a  function 
of  the  recirculation  degree.  It  is  worth  noting  that  the 
sum  of  the  two  efficiencies  approaches  unity  as  in 
the  type  of  burner  works  through  (mainly  irradiative) 
heat  losses  towards  the  surroundings.  It  can  be  also 
observed  that  the  radiation  efficiency  increases  with 
increasing  recirculation  degree;  thus  the  burner  en¬ 
sures  better  performances  in  the  MILD  combustion 
regime  than  in  the  flame  mode. 


4.6.  Comparison  of  axisymmetric  and  3 -dimensional 
simulations 


As  mentioned  earlier,  simulations  of  recuperative 
MILD  combustion  burners  available  in  the  literature 


are  2-dimensional  simulations  [19,20].  Therefore  the 
presence  of  the  three  recirculation  windows  is  usually 
neglected. 

In  the  present  work  axisymmetric  simulations 
were  also  performed  and  compared  to  the  3 -dimensio¬ 
nal  simulations  in  order  to  estimate  the  effect  of  the 
recirculation  windows,  as  well  as  errors  arising  when 
they  are  not  taken  into  account. 

Axisymmetric  simulations  will  be  denoted  as 
“2D”  for  brevity,  even  if,  strictly  speaking,  they  are 
not  2-dimensional  simulations,  the  computational 
domain  corresponding  to  a  5°  angular  sector  (see 
description  of  computational  domain  given  in  Sec¬ 
tion  3). 

Fig.  11  indicates  that  recirculation  degrees  pre¬ 
dicted  through  the  axisymmetric  model  are  larger  than 
those  predicted  through  the  3 -dimensional  model  by 
15-20%  for  the  smallest  air  inlet  cross-sectional  ar¬ 
eas.  This  is  reasonable,  as  in  axisymmetric  models  the 
exhaust  gases’  recirculation  into  the  reaction  zone  is 
facilitated  by  the  larger,  i.e.,  360°  opening.  For  the 
largest  air  inlet  cross-sectional  areas,  corresponding 
to  the  flame  mode  operations,  there  is  little  difference 
between  the  3 -dimensional  simulation  and  the  ax¬ 
isymmetric  model;  this  may  be  partly  imputed  to  the 
small  value  of  predicted  recirculation  degrees,  which 
makes  differences  negligible. 

Because  of  the  larger  predicted  recirculation  de¬ 
grees,  2-dimensional  models  are  expected  to  underes¬ 
timate  NO  emissions,  as  confirmed  by  Fig.  12. 

However,  results  obtained  with  the  two  approaches 
(2D  vs  3D)  could  be  fitted  by  the  same  regression 
curve  when  reported  as  functions  of  the  recircula¬ 
tion  degree  (see  Fig.  13),  suggesting  that  the  relation 
between  NO  formation  and  exhaust  gas  recirculation 
remains  the  same. 

The  best  regression  law  was  found  to  be  the  Boltz¬ 
mann  sigmoidal  regression 
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Fig.  12.  NO  emissions  as  a  function  of  the  air  inlet  cross- 
sectional  area  predicted  with  the  axisymmetric  and  3 -dimen¬ 
sional  models.  Burner  load  Q[n  =  10.42  kW  (runs  7-12  of 
Table  1). 
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Fig.  13.  NO  emissions  as  a  function  of  the  recirculation 
degree  predicted  with  the  axisymmetric  and  3 -dimensional 
models.  Burner  load  Qm  =  10.42  kW  (runs  7-12  of  Ta¬ 
ble  1). 
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Fig.  14.  Radial  profiles  of  temperature  for  different  recircu¬ 
lation  degrees  and  for  axisymmetric  and  3 -dimensional  mod¬ 
els  at  different  axial  coordinates:  (a)  x  =  150;  (b)  250  mm. 
Burner  load  Q{n  =  10.42  kW  (runs  7,  10,  11,  and  12  of  Ta¬ 
ble  1). 

5.  Conclusions 
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^NO_l  ~  ^NO_2 
l  _|_  6^R-^Rq/ Ak 


+  ^NO_2: 


(12) 


with  i  =  0.00553,  2  —  0.00001,  — 

—22.3%,  and  Ak  =  29.8%,  and  a  root  mean  squared 
error  approximating  unity. 

Therefore  axisymmetric  simulations  could  be  per¬ 
formed,  allowing  lower  computational  cost,  and  pre¬ 
dicted  NO  corrected  by  assuming  15-20%  larger  re¬ 
circulation  degrees. 

Nevertheless,  it  is  worth  noting  that  a  more  accu¬ 
rate  observation  of  the  burner  temperature  and  ve¬ 
locity  distribution  indicates  that  axisymmetric  and 
3 -dimensional  simulations  differ  substantially.  This 
is  evident  from  Fig.  14,  which  compares  radial 
temperature  profiles  obtained  with  the  two  models 
for  different  axial  coordinates  and  recirculation  de¬ 
grees. 


A  numerical  investigation  through  computational 
fluid  dynamics  of  a  recuperative  MILD  combustion 
burner  operating  in  MILD  combustion  mode  has  been 
presented. 

The  numerical  model  contained  two  subroutines: 
the  first  subroutine  models  the  air  preheater,  whereas 
the  second  subroutine  is  aimed  at  modeling  irradiative 
heat  losses  from  the  burner  toward  the  surroundings. 
Importantly,  in  this  manner  the  burner  model  was  self- 
sufficient,  which  means  it  did  not  need  any  assump¬ 
tions  (e.g.,  on  the  air  inlet  temperature). 

Because  of  its  complexity,  the  CFD  model  was  as¬ 
sessed  through  comparison  with  experimental  data, 
concerning  temperature  profiles  along  the  outer  sur¬ 
faces  of  the  burner  and  species  concentration  in  the 
flue  gases.  It  is  noteworthy  that  only  global  data 
and  measurements  taken  on  the  outer  surfaces  of  the 
burner  were  available,  as  the  burner  is  a  real  indus¬ 
trial  burner.  Local  data  would  be  desirable,  but  the 
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burner  does  not  allow  any  access  to  its  interior  for  the 
measurements.  Besides,  the  complexity  of  the  CFD 
model  and  the  need  for  subroutines  arose  partly  from 
the  wish  to  model  the  burner  without  any  additional 
assumptions,  trying  to  exploit  all  the  available  mea¬ 
surements  (e.g.,  temperature  along  the  radiant  tube). 

The  burner  was  observed  to  produce  about  30 
ppmv  of  NO  in  flameless  mode.  For  the  same  air¬ 
preheating  the  burner  would  have  produced  about 
1000  ppmv  of  NO  in  flame  mode.  In  addition,  it  was 
found  that  the  burner  performance  in  terms  of  thermal 
efficiency  improved  when  it  was  operated  in  MILD 
combustion  mode.  The  temperatures  on  the  outer  sur¬ 
face  of  the  burner  were  also  more  homogeneous,  en¬ 
suring  more  uniform  conditions,  which  are  beneficial 
for  many  process  such  as  steel  annealing  processes. 

A  further  reduction  of  the  air  injection  nozzles  was 
also  investigated  in  order  to  promote  recirculation  of 
exhaust  gases  by  increasing  the  air  inlet  jet  momen¬ 
tum  and  suppressing  further  NO  emissions.  In  partic¬ 
ular,  NO  emissions  of  just  10  ppmv  were  produced  by 
halving  the  air  inlet  cross-sectional  area. 

The  analysis  of  turbulence  and  chemical  timesca¬ 
les  confirmed  that  in  MILD  combustion  the  high-tem¬ 
perature  region  is  characterized  by  Damkohler  num¬ 
bers  lower  than  those  observed  in  the  flame  front  of 
the  burner  operating  in  flame  mode.  Therefore  the 
importance  of  the  turbulence-chemistry  interaction 
treatment  was  confirmed. 

Finally,  two  geometrical  models  of  the  burner  have 
been  compared:  a  3 -dimensional  model  and  an  ax- 
isymmetric  model.  This  latter  model  is  usually  pre¬ 
ferred  to  the  former  in  the  available  literature  for 
the  modeling  of  MILD  combustion  burners,  because 
of  its  lower  computational  cost.  It  was  found  that 
the  axisymmetric  model  overestimates  the  recircula¬ 
tion  degree  in  the  burner  by  15-20%,  and  thus  un¬ 
derpredicts  NO  emissions.  However,  for  engineering 
purposes,  a  satisfactory  prediction  of  NO  emissions 
could  be  achieved  by  extrapolating  data  from  the  NO 
concentration  versus  the  recirculation  degree  curve 
and  considering  a  lower  recirculation  degree  than  that 
predicted  through  the  axisymmetric  model.  However, 
temperature  and  flow  fields  obtained  with  the  axisym¬ 
metric  model  should  be  handled  carefully. 

In  conclusion,  the  careful  application  of  commer¬ 
cial  codes  with  detailed  analysis  of  different  model¬ 
ing  hypotheses  is  able  to  give  a  satisfactory  view  of 
the  most  important  effects  of  MILD  combustion.  In 
particular,  in  industrial  conditions,  where  the  burner 
structure  and  the  characteristics  of  the  experimen¬ 
tal  tests  do  not  make  it  possible  to  take  detailed  lo¬ 
cal  measurement  within  the  reacting  system,  the  ap¬ 
proach  followed  describes  the  main  features  of  the 
reacting  system,  at  least  from  an  engineering  point  of 
view. 
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